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Abstract. We study a certain random growth model in two dimensions closely 
related to the one-dimensional totally asymmetric exclusion process. The results 
show that the shape fluctuations, appropriately scaled, converges in distribution 
to the Tracy- Widom largest eigenvalue distribution for the Gaussian Unitary En- 
semble (GUE). 
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1. Introduction and results 

The shape and height fluctuations in many 2-d random growth models are 
expected to be of order N x , with \ = 1/3, if the mean of the linear size of the 
shape or the height is of order N. See [KS] for a review and [NP] for rigorous 
bounds on \ in first-passage percolation. 

In this paper we will consider a specific model. It can be given several proba- 
bilistic interpretations, as a randomly growing Young diagram, a totally asymmet- 
ric one dimensional exclusion process, a certain zero-temperature directed polymer 
in a random environment or as a kind of first-passage site percolation model. The 
model has the advantage that we can prove that % = 1/3 and also compute the 
asymptotic distribution of the appropriately rescaled random variable. Interest- 
ingly, the limit distribution that occurs is the same as that of the scaled largest 
eigenvalue of an N x N random matrix from the Gaussian Unitary Ensemble 
(GUE) in the limit N — > oo. The model in this paper has many similarities with 
the problem of the distribution of the length of the longest increasing subsequence 
in a random permutation where the same limiting distribution and % = 1/3 was 
found in [BDJ]. 

To define the model let w(i,j), G Z+, be independent geometrically 

distributed random variables, 

¥[w(i,j) = k] = (1 -q)q k , keN, 

where < q < 1. Let Hm,n be the set of all up/right paths 7r in I? + from (1, 1) 
to (M, N), i.e. sequences (ik,jk), k = 1, . . . , M + N — 1, of sites in Z^_ such that 
(iiji) = (1,1), (iM+N-ijM+N-i) = (M,N) and (i k+1 J k+1 ) - (i k Jk) = (1,0) 
or (0, 1). Define the random variable 

G(M,N)= max V w(i,j). (1.1) 

(l,j)67T 

We also define the closely related random variable 

G*(M,N)= max V w*(i,j), 

0,j)67T 

where w*(i,j) = w(i,j) + 1, so that P[w*(i,j) = k] = (1 - q)q k ~ 1 , k > 1. Clearly, 



G*(M,iV) = G(M, N) + M + N- 1, 



(1.2) 
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since all paths have the same length. Using this random variable we can define, 
for each t > 0, a random subset of the first quadrant by 

A(t) = {(M, N)eZ 2 + ; G*(M, N)<t} + [-1, 0] 2 . (1.3) 

From the definition of G*(M, N) and the fact that we consider up/right paths it 
follows that A(t) has the form 

Ul^[k-l,k} x [0,A fc ] 

for some integers Ai > A2 > . . . > A r > 1, so we can think of A(t) as a Young 
diagram A = (Ai, . . . , A r ). If we think of t e N as a discrete time variable, A(t) is 
a randomly growing Young diagram. Let d*A(t) be those unit cubes adjacent to 
A(t) that can be added to A(t) so that it is still a Young diagram, i.e. each cube in 
d*A(t) must have a cube in A(t) or R 2 \ [0, oo) 2 immediately below and to the left 
of it. The fact that the w*(i,j) :s are independent and geometrically distributed 
random variables implies that A(t + 1) is obtained by picking each cube in d*A(t) 
independently with probability p = 1 — q and adding those cubes that were picked 
to A(t). (Recall that V[w*(i,j) = k+l\w*(i,j) > k] = F[w(i,j) =/],/> 0, the lack 
of memory property.) The starting configuration is A(0) = and d*A(0) = [0, l] 2 . 
In this model G*(M, N) = k means that the box [M - 1, M] x [N - 1, N] is added 
at time k. This growth model has been considered in [JPS]. 

This randomly growing Young diagram can also, equivalently, be thought 
of as a certain totally asymmetric exclusion process with discrete time, compare 
[Ro] or [Li], p. 412. Let C(t) = <9([0, oo) 2 \ A(t)) and note that C(t) consists of 
vertical and horizontal line segments of length 1. To each vertical line segment we 
associate a 1 and to each horizontal line segment a 0. If we read the numbers along 
C(t), starting at infinity along the y-axis and ending at infinity along the x-axis, 
we get an infinite sequence X(t) = (. . . , x-\(t), xo(t), xi(0), £2(0), . . .) of 0's and 
l's, starting with infinitely many l's and ending with infinitely many 0's; we let 
xo be the last number we have before passing through the line x = y. We can 
think of X(t) as a configuration of particles, where Xk = 1 means that there is a 
particle at k, whereas Xk = means that there is no particle at k. The stochastic 
growth of A(t) described above corresponds to the following stochasic dynamics 
of the particle system. At time t each particle independently moves to the right 
with probability 1 — q provided there is no particle immediately to the right of 
it. Otherwise it does not move. The starting configuration is Xk(0) = l(_ 00j o](fc). 
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In this particle model G*(M, N) = k means that the particle initially at position 
— (N — 1) has moved M steps at time k. 

Our first result concerns the mean and large deviation properties of G(M, N). 

Theorem 1.1. For each q e (0, 1) and 7 > 1, 

lim i-E[G([ 7 iV] , AT)] = (1 1 )2 - 1 = o;( 7 , g). (1.4) 

N^oo iv 1 — (/ 

Also, Cr([7iV], A/") has tie following large deviation properties. There are functions 
i(e) and £(e) (which depend on q and 7), so that, for any e > 0, 

hm i^logP[G([7iV],iV)<iV(a;(7,g)-6)] = -£(e) (1.5) 

N^oo iv 

and 

hm llogP[G([ 7 iV],A0 > JV(w(7,g) + e)] = -i(e). (1.6) 

N^oo iv 

The functions £(x) and i(x) are > if x > 0. 

Note that the existence of the limit (1.4) follows by a subadditivity argument, 
so it is the explicit form of the constant that is interesting. The large deviation 
result (1.6) has been obtained in [Se2]. The theorem will be proved in section 2. 

The theorem implies that jA(t) has an asymptotic shape A as t — > 00, in 
the sense that given any e > 

(l-e)Ao C -A{t) C (l + e)A 

for all sufficiently large t. It follows from the definition of A(t), (1.3), and theorem 
1.1 that 

A = {{x, y) E [0, oo) 2 ; y + 2^xlj + x<l-q}. 

The boundary of Aq consists of two line segments from the origin to (1 — q, 0) and 
(0, 1 — q) and part of an ellipse that is tangent to the x- and y-axes. 

We now want to understand the fluctuations of A(t) around its asymptotic 
shape Aq, i.e. the fluctuations of G([yN], N) around Nuj^^q). Before we can 
formulate the result we need some preliminaries. Let Ai (x) be the Airy function 
defined by 

Ai(x) = — / e m+isf/3 + ix(t+is) dt 
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where s > is arbitrary. Consider the Airy kernel 

Ai W Ai-(»)-Ai' W Ai(») 
x-y 

as an integral kernel on L 2 [s, oo). The Fredholm determinant 

F(s) = det(I-A)\ L . [Sj00) =J2 L -n L / det(A(^, Xj )% =1 d k x (1.8) 

is a distribution function. It is the distribution function of the appropriately 
scaled largest eigenvalue of an N x N random matrix from the Gaussian Unitary 
Ensemble (GUE) in the limit N — > 00, the Tracy- Widom distribution, see [TW1]. 
The distribution function F(s) can also be defined using a certain Painleve II 
function, 

/■oo 

F(s) = exp[-/ (x - s)u(x) 2 dx], (1.9) 

J s 

where u(x) is the unique solution of the Painleve II equation 

u" = 2u s + xu, 

with the asymptotics u(x) ~ Ai (x) as x — > 00. The fact that the expressions (1.8) 
and (1.9) are equal is proved in [TW1]. 

Theorem 1.2. For each q e (0, 1), 7 > 1 and s£l, 



where 



^(7, q) = ^l 1 ^ ^ + v^) 2/3 (l + v^7) 2/3 - (1-11) 



The theorem will be proved in section 3. We have not proved convergence 
of the moments of the rescaled random variable, see remark 2.5. This theorem 
should be compared with the result obtained in [BDJ], that if In(&) is the length 
of a longest increasing subsequence in a random permutation a G Sn (all AH 
permutations have the same probability), then 

lim F[(VN)- 1/3 (£ N (a) - 2VN) < s] = F(s). (1.12) 

Note that in both cases we have the same exponent 1/3, the standard deviation is 
~ (mean) 1 / 3 

The proofs of theorems 1.1 and 1.2 are based on the following result which 
will be proved in section 2. 
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Proposition 1.3. For any M > N > 1, 

w^=^ e n o.-^n( Ai+ r JV )^' 

M ' iV heN N \<i<j<N i=l V 4 7 

max{/ii}<t+Af-l 

(1.13) 

where Zm,n is the normalization constant (partition function). 

This remarkable formula should be compared with the formula for the distri- 
bution function for the largest eigenvalue, A max , of an iV x N random matrix from 
GUE, 

1 r - N 

P[A max <*] = —/ [T (x t - x 3 f TT e~ 2Nx U N x. (1.14) 

There is a clear similarity between the two expressions, so we can use the ideas 
developed to investigate (1.14). Just as the right hand side of (1.14) can be written 
as a Fredholm determinant, so can the right hand side of (1.13). The kernel for 
(1.13) is the Meixner kernel, 

Km,n(x,v) 

= KN _, M^xW^jy) - M^WMsjy) (<(a;)<(y)) i/2 > (L15) 
kn x — y 

where Mn(x) = knx n + ... are the normalized orthogonal polynomials with 
respect to the discrete weight, K = M — N + 1, 

w q K (x)= (^ + I ^~ 1 ^jq x , xeN. (1.16) 

This Meixner kernel also appears in the recent paper [BO]. The polynomial Mn{x) 
is a multiple of the classical Meixner polynomials m^' q (x). Using the explicit gen- 
erating function for the Meixner polynomials, see [Ch], the appropriate asymp- 
totics of the kernel (1.15) can be analyzed. This will be done in section 5. 

Let u(i,j), G be independent exponentially distributed random 

variables with parameter 1. Let H(M, N) be the analogue of G(M, N) for these 
random variables, i.e. 



H(M, N) = max{ j) 5 * e U m ,n}- 



(1.17) 
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We can consider the related stochastically growing Young diagram and totally 
asymmetric exclusion process just as in the geometric case, where we now have 
continuous time. This simple exclusion process is exactly the one considered by 
Rost, [Ro], see also [Li]. In this process X(t) = (VkiWkL-oc e {A 1} Z the initial 
configuration is l(_ 00j o](fc) and a particle (rjk = 1) jumps with exponential rate to 
the right one step provided there is no particle at k + 1 (rjk+i = 0). By taking the 
q — > 1 limit in (1.13) we obtain 

Proposition 1.4. For any M > N > 1, t > 0, 

i r N 

¥[H(M, N)<t] = / TT ( Xi - x 3 f T\ xf- N e~^d N x. (1.18) 

Z M,N J[0,t]» 1 <i < j< N fA 



PROOF: If Xl is geometrically distributed with parameter 1 — 1/L, then L~ 1 Xl 
converges in distribution to an exponential random variable with parameter 1. 
Since G(M, N) is a continuous function of the w(i,j) :s, proposition 1.3 gives 



F[H(M, N) < t] 

L^oo Z M ,N i ^ ~7 \ h i J 



(*) l<i<j<7V i=l 

rJV 2 7 , M — N 1 .1 

lim £ _y FT ( ^Z^)2TT e -^+°(i) TT (^±^) 

l^oo Z M ,7v M - iV ! .11 1 L 11 11 L 

' v (*) 1<i<j<N i=l k=l 



tf / n (^-^) 2 n 

'M,iV -/[O.tp 1<i<j<N ? = i 



iV 



where (*) means summation over all h G N N such that max{/ij} < [Lt] + N — 1. 

Remark 1.5. The right hand side in (1.18) is the probability that the largest 
eigenvalue in the Laguerre ensemble is < t. It occurs in the following way. Let A 
be an N x M rectangular matrix (N < M) with entries that are complex Gaussian 
random variables with mean zero and variance 1/2. Then the right hand side in 
(1.18) is the distribution function for the largest eigenvalue of AA* , see [Ja]. 

Theorem 1.6. For each 7 > 1, 

lim ±-E[H(tyN], N)] = (1 + ^) 2 , (1.19) 
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and there are functions z*(e) and £*(e) (which depend on 7), so that for any e > 0, 
lim -Ll og F[H([ 1 NlN)<N((l + ^) 2 -e)] = -L(e) (1.20) 

TV— »oo iv 



and 



lim llogP[ J ff([ 7 AT],iV)>iV((l + V7) 2 + e)] = -^(e). (1.21) 



Furthermore, assume that = 0(N 1 ^ 3 ) as N — > oo and pick djv so that rf^ — 
(1 + l/y/7)a N = oiN 1 / 3 ) as N -> oo. Then, for each 7 > 1, 

iV^oo L 7-V6(l _|_ ^4/3^1/3 - J V 7 V ; 

Proof: For the proof of (1.19) to (1.21) see remark 2.3. Write c = (1 + v /t") 2 and 
p = 7 _1 / 6 (1 + y/^f) 4 ^ 3 - Then, by proposition 1.4, 

F[H(iN + a N , N) < cN + d N + piV 1/3 s] 



1 /■ N 
/ A{x) 2 T\xJ N e- x ^d N x, 



Z'~(N+a N ,N J[0,cN+d N +pN 1 / :i s} 

where A(x) = rii<i<j<iv( a '.? ~ x i) an d gin = (7 — l)iV + aAr. By a standard argu- 
ment, see [Me], ch. 5, [TW3] or section 3, this equals the Fredholm determinant 



J2 ( -TT- I detipN^K^icN + dN + pN^cN + dN + pN 1 / 3 ^^^ 
k=o k - J i s ^) k 

(1.23) 

where 



is the Laguerre kernel. Here, 

/ n' \ 1/2 

C(*) = ( (-l)^n(^) = + • • • 

are the normalized associated Laguerre polynomials, 

\ x )^m (x)x ^ cfe = 5 nm . 
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From asymptotic formulas for these polynomials it follows that 

lim KZ N (cN + d N + pN 1 /%cN + d N + pN 1 / 3 ^) = A(£,r]). (1.24) 

N^oo 

This can be proved in the same way as the corresponding results for Meixner 
polynomials, see sections 3 and 4, by using the integral representation 

e x f e~ xz z n + a 



L " {x) -2^J c (z-ir+i dz > 

where C is a circle surrounding z = 1. Using (1.23), (1.24) and some estimates 
(compare lemma 3.1) we obtain 



TV— >oo 



lim P[#( 7 iV + a N , N) < cN + d N + piV 1/3 s] 
E ij TT I det(A(6,^))t=i^ = F{s). 



k=0 



We will not present all the details since they are similar to the proof of theorem 
1.2. 

Using this result we can get a fluctuation theorem for Rost's totally asym- 
metric simple exclusion process defined above. The random variable H (N, M) is 
the first time at which the particle starting at —(N — 1) has moved exactly M 
steps to the right. If we define Y(k, t) = J2j>k Vj(t) to be the number of particles 
to the right of k at time t. Then Y{k,t) > m means that the particle that starts 
at — m has moved > m + k + 1 steps at time t. Hence 

F[Y(k, t)<m] = l- F[H(m + k + 1, m + 1) < £]. 

Using this relation and (1.19) to (1.21) we obtain the follwing result first proved 
by Rost, [Ro], 

almost surely as t — > oo, \u\ < 1. Now, using (1.22) it is fairly straightforward to 
show the following result. 

Corollary 1.7. For each u e [0, 1), 

hm P[Y M , t) <l(l-uf + |^P^ 1/3 ] = 1 " F(-0- 
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Remark 1.8. We can interpret theorems 1.1 and 1.2 (and analogously theorem 
1.6) as a result for a kind of zero-temperature directed polymer or equivalently a 
directed first-passage site percolation model in the following way. 

Let Sk be the simple random walk in Z starting at at time and ending at 
at time 2N + 2. Denote the set of all possible paths by Vn- Let v(i,j), G Z 2 
be independent, identically distributed random variables, and let (3 > 0. On Vn 
we put the random path probability measure 

2N 

Q ( N [S] = - F eM-Pj2 v ^ S ^ 

N k=l 

S E Vn, where C^r is the normalization constant. This measure describes a di- 
rected polymer (S) fixed at both endpoints at inverse temperature (3 in the random 
environment given by the v(i,j) :s, see [Pi]. The free energy is — logC^, and 
in the zero temperature limit (3 — > oo this becomes 

27V 

ZeVN k=i 

the ground state energy. By rotating the coordinate system by the angle — it/ 4 
it is seen that (1.25) can, equivalently, be thought of as a first-passage time in a 
directed first passage site percolation model. Let u{i,j), (z,j) G R+, be indepen- 
dent, identically distributed random variables (with the same distribution as the 
v(i,j) :s). Then F^ s has the same distribution as F(N, N), where 

F(M,N)= min V] u(i,j). 

0,j)67T 

(The u(i,j) :s are usually thought of as passage times and F(M, N) is the minimal 
flow time from (1,1) to (M, N). Hence it is natural to assume that u(i,j) > 0, 
but this will not be the case below.) We can define a random shape 

B(t) = {(M, N)eZ 2 + ; F(M, N) < t} + [-1, 0] 2 . 

Set u(i,j) = a — w(i,j), where a > a min = (1 — q)~ 1 (q + y/q) (this condition on 
a ensures that B{t) will grow); w(i,j) are the geometrically distributed random 
variables considered above. Then clearly, 



F(M, N) = a(M + N - 1) - G(M, N), 



(1.26) 
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Let 7 > 1, set £(7) = (1 + 7 2 )~ 1/2 (7, 1), a unit vector and [nx(^)\ = ([JV7], JV), 
([•] the integer part, where N = [(1 + 7 2 ) _1 / 2 n], so that [^£(7)] is a lattice site 
near nx(^). Let T n ( 7 ) be the first time s > for which B(s) reaches [71^(7)], 

T n ( 7 ) = inf{s > 0; [nx( 7 )] G 

Clearly, by the definition of -B(s) and equation (1.26), 

T n ( 7 ) = a([>yN] + N - 1) - G([ 7 iV], iV), 

where AT = [(1 + 1 2 )- 1 / 2 n]. 

Theorem 1.1 implies that for each q e (0, 1) and 7 > 1, 

hm iE[T n ( 7 )] = _L=[a( 7 + 1) - (1 ± ^ + 1] = /*( 7 ). 
n^oo n Ji -|- 7 2 1 — Q 

Also, T n ( 7 ) has large deviation properties similar to those for G([yN],N). Using 
this result we can compute the asymptotic shape of B(t). It follows from theorem 
1.2 that 



L (l + 7 2)-l/6 p ( ?)7 ) n l/3 

as n — > 00. 



<s]-l-F(-s), 



Conjecture 1.9. Is the result for G([yN], N) limited to geometric and expo- 
nential random variables? Normally, we expect limit laws for appropriately scaled 
random variables to be independent of the details. It is therefore natural to 
conjecture that if the w(i,j):s are i. i. d. random variables with some suit- 
able asumptions on their distribution, then there are constants a and b so that 
(G([yN], N) — aN)/bN 1 / 3 converges to a random variable with distribution F(s). 
By remark 1.8 this leads to a related conjecture for directed first-passage site 
percolation. 

2. The Coulomb gas 



2.1 Combinatorics. 

The key combinatorial ingredient is the Knuth correspondence introduced 
in [Kn]. It generalizes the Schensted correspondence [Sc] which is used in [BDJ]. 
Write [N] = {1, . . . , N}. Let Mm,n denote the set of all Mx N matrices A = (a^-) 
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with non-negative integer elements, and let M k M N be the subset of those matrices 



that satisfy YlfiLi ^27=i a ij = k. A two-rowed array 



^N 



^l 
ji 



jk 



is called a generalized permutation if the columns Q r ) are lexicographically or- 
dered, i.e. either i r < i r+ i or i r = i r+ i, j r < j r +i- There is a one-to-one 
correspondence between the set S 1 ^ N of all generalized permutations of length k, 
where the elements in the upper row come from [M] and the elements in the lower 
row from [AT], and M k M N defined by a — > f(a) = A = (ay), where 



dij = #times 



occurs m a. 



We say that ( l . ri ),..., f Vm ), ri < r 2 < . . . < r m is an increasing subsequence 
in cr if < j 2 < . . . < j rm - Let £(cr) denote the length of a longest increasing 
subsequence in a. 

Example. The generalized permutatation 



1 1 1 2 2 2 3 3 4 4 
1 2 2 2 2 2 1 2 1 3 



corresponds to 



(1 2 OX 

3 

1 1 
Vl 1/ 

A longest increasing subsequence is 12222223 so £(a) = 8. 

Recall from section 1 that Hm,n denotes the set of all up/right paths 7r from 
(1, 1) to (M, N) through the sites with 1 < i < M, 1 < j < N. 



Lemma 2.1. For each A e M k M N , 



max{ ^ Oy ; n G TI m ,n} = i(f 1 (A)). 



(2.1) 



PROOF: This is clear from the definitions. That we go to the right corresponds 
to the fact that i ri < . . . < i Tm and that we go up corresponds to j ri < . . . < j Tm 
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(the upper row gives row indices whereas the lower row gives column indices in 
the matrix). 

Now, Knuth has defined a one-to-one mapping from the set N to pairs 
(P, Q) of semi-standard Young tableaux of the same shape A, which is a partition 
of k, A h k, where P has elements in [N] and Q has elements in [M]. (More 
information on Young tableaux can be found in [Sa] and [Fu].) This correspondence 
has the property that if a — > (P, Q) and P, Q have shape A, then £(a) = the length 
of the first row, Ai, in A. Consider G(M, N) defined by (1.1). The M x N matrix 
W = (w(i,j)) is a random element in M.m,n- Let 

M N 

S(M,N) = J2Y, w ^i) 
i=i j=i 

and 

PM, N (t)=nG(M, N)<t]. 

Then, 

oo 

PM, N (t) = ^F[G(M,N) < t\S(M,N) = k]F[S(M,N) = k}. (2.2) 

fc=0 

For a fixed A E M k M N we have 

n{A}] = - q)q a - = (1 - q) MN q\ 

since J2i j a ij = k- We have proved 

Lemma 2.2. The conditional probability F[-\S(M, N) = k] is the uniform distri- 
bution on M\j N . 

This lemma is the reason that we choose the w(i,j) :s to be independent and 
geometrically distributed. Note that 

F[S(M, N) = k} = #M k M , N {l - q) MN q k . (2.3) 

Let L(A, M, N) denote the number of pairs (P, Q) of semi-standard Young tableaux 
of shape A, such that P has elements in [N] and Q has elements in [M]. Combining 
lemma 2.1, lemma 2.2 and the Knuth correspondence we see that 

F[G(M, N) < t\S(M, N) = k] = \ £ L(X,M,N). (2.4) 

W JVl M,N Ahfc)Al < t 

To compute L(A, M, N) we use 



15 



Lemma 2.3. The number of semi-standard tableaux of shape X and elements in 
[N] equals 

At -Xj+j-i 



n 



l<i<j<7V J 

PROOF: We have two formulas for the Schur polynomial in N variables associated 
with the partition A, [Sa], [Fu], 

^-^ rp det(x^ l+N l )i<i,j<N 

S\(X) — y X = — Tf : , 

1 } V detfcif-)!^* 

where the sum is over all semi-standard A-tableaux T with elements in [N] and 
N with rrij equal to the number of times j occurs in T. Hence, 
evaluating the Vandermonde determinants, 

ry.Xi-Xj+j-1 _ 1 

l<i<j<N 

where r = X)il=i(* ~~ 1)<V The number of semi-standard tableaux with elements 
in [N] equals 

Xi -Xj+j-i 



s x (l,l,...,l)=\ims x (l,x,...,x N - 1 )= 11 



x—>l 1 — I 

l<i<j<N J 

This completes the proof of the lemma. 

It follows from lemma 2.3 that 

L(X,M,N)= TT h^i±A^i TT Xi-Xj+j-i^ (25) 
j-j. j — i - LJ - j — i 

l<i<j<M J l<i<j<N J 

We assume from now on that M > N, the other case is analogous by symme- 
try. Since the numbers in the columns in P and Q are strictly increasing we must 
have Xi = if N < i < M. Hence 

m, n) — n ( A --^-' )n n (^)- 

l<i<j<M V J 7 i=l j=N+l V J 7 

Let /ij = Aj + N — j, j = 1, . . . , N, so that hi = Ai + N — 1, /ijv = A^v > and 
foi > /i2 > • • • > ^at- Then 

l(x,m,n )= n %e^n n 

l<j<j<7V VJ ; i=l j=JV+l J 
j=0 J V Jy l<i<j<N i=l 1 
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The condition J2f=i -\? = ^ translates into Yl^Li hj = k + N(N — l)/2 and Ai < t 
to h\ < t + N — 1. By (2.2), (2.3) and (2.4) we have 

oo 

PM,iv(t) = - ?) M V E £(A,M,JV) 

fc=0 Ahfc,Ai<fc 

and inserting (2.6) yields 

N-l 



( i n \MN iv -- 
m l 1- ^ -N(N-l)/2 TT 



AT! A± j!(M-iV + j)! 



k=0 heK N l<i<j<N i=l l ' 

hi=k+N(N-l)/2 
max{hi}<t+N-l 

where we have used the symmetry under permutation of the hi :s. Summing over 
k gives all the possible values oi^hi, so we obtain 



1 N 
PM, N (t) = - H (hi-hj) 2 ]l w M-N + i(hi)- (2-7) 

J M,N 



hen N i<i<j<N i=i 

max{hi}<t+N-l 



where w q K {x) is given by (1.16) and 

JV-l 

Z M ,N = q^-^il - q)- MN H j\(M - N + j)l (2.8) 

3=0 

This proves proposition 1.3. 

2.2 The large deviation estimate. 

In order to investigate the location of the rightmost charge in (2.7) and prove 
large deviation formulas we rescale the discrete Coulomb gas (2.7). Let M = pyiV], 
7 > 1 fixed, and K = K(N) = [ 7 iV] - M + 1. Set = ^N, A N (s) = {x E 
An; x < s} and 

V™(t) = -±logw q K{N) (Nt), t>0. 
Using Stirling's formula we see that 

lim V$ q (t) =tlog--(t+7-l)log(t+7-l)+tlogt + (7-l)log(7-l) = V^' q (t) 
N^oo q 

(2.9) 
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uniformly on compact subsets of [0, oo). (We will often omit the superscripts 7 
and q.) Rescaling the variables in (2.7) by setting hi = Nxi, Xi E An we see that 
(2.7) can be written 

PN(t) = PM(N)Mt) = Zn ^ + 1 ~^\ (2.10) 

where 

JV 

Z N (s)= A 7V (x) 2 exp(-iV^y Ar (x J )) (2.11) 

xeA N (s) N 3=1 

and Zn = Zn(oo). Here An(x) = rii<i<j<jv( a '.7 ~~ x i) * s the Vandermonde 
determinant. 

When investigating the large deviation properties of pjv(£) we may just as 
well consider more general confining potentials Vn- Assume that Vn '■ [0, 00) — > R, 
N > 1, satisfy 

(i) Vn is continuous, N > 1. 

(ii) There are constants £ > 0, T > and N > 1 such that 

V N (t) > (l+0log(t 2 + l) (2.12) 

for t > T and N > N . 

(iii) Vjv(^) - ► uniformly on compact subsets of [0, 00). 
Set for x EAfJ and /3 > 

Qm,n{x) = \ A M (x)f JJexp( — —V N (xj)). 

i=i 

(This M is not the same as the previous M.) Define the partition functions 

ZM,N{t) = ^ Qm,n(x), 

xeA N (t) M 

Zm,n = Zn,m{oo) and the probability measure 

Pm,n[B] = — *— V Qm,n(x), 

' x£B 

B C N M . We are interested in the distribution of the position of the rightmost 
charge, maxi<fc<M Xk- Its distribution function is given by 

F M , N (t) = P M ,7v[maxa: fc <t}= Z f N ^\ (2.13) 
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(If M = N we write F N (t).) 

In order to formulate the large deviation results for Fjv(t) we need some 
results from weighted potential theory, [ST]. The results we need differ from the 
usual ones since we are interested in the continuum limit of a discrete Coulomb 
gas, so that the particle density of the rescaled gas is always < 1. Hence, the 
equilibrium measures will be absolutely continuous with a density <fi satisfying 
< < 1. Let A s denote the set of all <fi G -^ 1 [0, s) such that < < 1 and 
J Q S = 1, 1 < s < oo. Given V : [0, oo) — > R, continuous and such that there is a 
5 > and a T > such that 

> (1 + 5) log(t 2 + 1) (2.14) 

for t > T, we set 

*v(s, y) = log | a; - y\~ l + \v(x) + \v(y) 



and 



I v[4>]= / / k v (x,y)(f)(x)(f)(y)dxdy, 
Jo Jo 



for G »4 S . 

The proof of the next proposition is similar to the corresponding result in 
weighted potential theory. See [DS] and also [LL] where a very similar problem is 
treated. 

Proposition 2.1. For each s G [1, oo] there is a unique 4> s v G A s such that 

inf IvM=Iv[#r] = F£. 

Tie extremal function (j)y has compact support. (If s = oo we will drop the 
superscript. ) 

Let 6y = sup(supp0\/) be the right endpoint of the support of 4>v- Set 
J{t) = for t < b v and 

/•OO 

J(t) = M k v (r,x)(p v (x)dx- F v (2.15) 

for t>by. Also, set 

L(t) = i(i^-Fv0 

for t > 1. The next theorem gives the large deviations for the distribution function 
F N (t) defined by (2.13) 
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Theorem 2.2. Assume that V^it) satisfies the assumptions (i) - (Hi) above. Then 

jfco^a logF N (t) = -0L(t) (2.16) 

for anyt > 1 and L(t) >0ift<by- Assume furthermore that J(t) > fort > by- 
Then 

lim ±log(l -F N (t)) = -(3J(t) (2.17) 

for all t. 

We postpone the proof to section 4. 

Remark 2.3. The same result is true for a continuous Coulomb gas on R with 
density 

-^\A N (x)feM-^J2nxj)), (2-18) 

Z N j=l 

on M. N , which occur in random matrix theory. The choice (3 = 2 and V(t) = 
2t 2 corresponds to the Gaussian Unitary Ensemble (GUE), compare (1.14). We 
assume that V is continuous and satisfies (2.14). In this case A s is replaced by 
M.i(s), the set of all probability measures on (— oo, s), and 4>v{x)dx is replaced by 
the equilibrium measure djiy (t) , see [Jo] . The proof is essentially the same. The 
formula (2.16) for certain V is a consequence of the result in [BG], see also [HP]. 
Also, (2.17) has been proved in the case V(t) = t 2 /2 in [BDG]. If we take (2.18) 
on [0, oo)^ with (3 = 2 and V(t) = —(M/N — 1) \ogt + t we get the measure in 
(1.18), and in this way we can prove (1.19) to (1.21). 

We can now apply theorem 2.2 to the model we are interested in. It is straight- 
forward to verify that V^' q satisfies the conditions (i) - (iii) with limiting external 
potential V' y,q (t). Write by^^ =6(7, q). The computation of <$>yi, <? will be outlined 
in section 6. We have 

(i + Vrf) 2 



6(7, q) = 



1-q 



If 7 > 1/q, then 



2 

<l>vy<i(t) = v(-(t — a) — 1), a <t < b, 
c 



where a = > c = ^(7 5 o) — a an d 

/ \ 1 r / Dx + 1 , . Bx + 1 . . 

v(x) = — arctanl - =— = — arctanl - — : ) , (2.19 
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B = (7 + q)/2y/rf, D = (1 + q^)/2^/qj. If 7 < V?, then, 
where 



1, if0<t<a 
u(§(*-a)-l), ifa<t<6, 



v(x) = — 7r — arctan ; — arctan ; 2.20 

with a, c, £?, D as before. 

We will not discuss the explicit form of the lower tail rate function. The upper 
tail rate function is given by 

with c, £?, D as above and x = 2{t — a)/c—l. Using this formula we can show that 
(see section 6) there are constants c\ > and C2 > so that 

r Cl 5 3 / 2 if < 5 < 1 

and 

J(b + 5)= — q > 7 S 3 / 2 + 0(5 5 / 2 ). (2.23) 

1 ' 3^/4(^+^(1 + ^) ^ J ^ ' 

In particular J(i) > if t > 6(7, q). 

From (2.10), (2.13) and theorem 2.2 we obtain 

lim -L log PJV (iVt) = -2L(t + 1) (2.24) 

N^oo iv 

and 

lim llog(l- Pjv (iVt)) = -2J(t+l) (2.25) 

JV— >oo iv 

for each £ > 0. These formulas imply theorem 1.1 with £(e) = 2L(by — e) and 
i(e) = 2J(6 y + e). By theorem 2.2 and (2.22) we have i(e) > and £(e) > if 
e > 0. 

By a superadditivity argument, the limit (2.25) actually gives a large deviation 
estimate for all N, compare [Sel]. 
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Corollary 2.4. For allt>0 and N > 1, 

1 -p N (Nt) < exp(-2NJ(t + 1)). (2.26) 

Proof: For 1 < M x < M 2 and 1 < JVi < A^" 2 we let G[(M 1 ,N 1 ),(M 2 ,N 2 )] 
denote the maximum of j) €7r over all up/right paths from (Mi, iVi) to 

(M 2 , iV 2 ). Note that if 1 < M x < M 2 and 1 < JVi < jV 2 , then 

(i) G[(Mi + l,A^i + l),(M 2 ,A^ 2 )] andG[(l,l),(M 2 -Mi,AT 2 -Ar 1 )] are identically 
distributed. 

(ii) G[(l, 1), (Mi, Nij\ and G[(M 1 + 1, AT X + 1), (M 2 , AT 2 )] are independent. 
Since [27 AT] > 2[yN], we have 

(iii) G[([yN] + 1,N+ 1), ([2 7 AT], 2AT)] > G[([yN] + 1,N+ 1), (2 [7 AT], 2 AT)]. 
Write OatW = 1 - p N (Nt) = F[G((1, 1), ([ 7 AT], iV)) > iVt]. Then, by (i) and 

(iii), 

aiv(t) < P[G(([ 7 A r ] + 1, AT + 1), ([2 7 AT], 2 AT)) > AT£] 

and hence, by (ii), a7v(t) 2 < a 2 Ar(t). Repeated use of this inequality yields 
A^ -1 logout) < (2 A: A r )~ 1 loga 2 fe7v(^) and by letting k — > 00 and using (2.25) 
we find N" 1 logaAr(t) < -2J(t + 1). 

Remark 2.5. We cannot prove convergence of the moments of the rescaled ran- 
dom variable in theorem 1.2 since we have no finite N estimate of P[G([7A r ], N) — 
ujN < — sN 1 / 3 ] for s > large. This would require an estimate of the finite N 
Fredholm determinant. In the other direction we can use the estimate in corollary 
2.4. The same remark applies to theorem 1.6. 

Remark 2.6. In [BR] it is proved by Baik and Rains that if we consider permuta- 
tions with certain restrictions we can get the Tracy-Widom distributions for GOE 
and GSE as limiting laws for longest increasing and decreasing subsequences. By 
considering a restricted geometry we can obtain the Tracy-Widom distribution for 
GOE, [TW2], also in the present setting. Let w(i,j) , 1 < % < j be independent 
geometrically distributed random variables, F[w(i,j) = k] = (l — q)q h for 1 < % < j 
and ¥[w(i,i) = k] = (1 — ^Jq)q k ^ 2 for % > 1. Set w(i,j) = w{j,i), if i > j > 1, 
so that A = (w(i,j)) is a symmetric matrix. The Knuth correspondence maps A 
to a pair of semistandard Young tableaux (P, Q) with Q = P, i.e. A maps to a 
single semistandard Young tableaux, see [Kn] or [Fu]. Let II^ 1 ^ be the set of all 
up/right paths from (1, 1) to (N, N) in {(2, j) G ; 1 < i < j}, i. e. in a triangle, 
and set 

F(A0=max{ £ w(i,j) ; tt e 11%™}. 

(i,j)€iv 
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Now, we also have 

F(N) = max{ ^ w(i,j) ; tt e n^w}, 

(i,j)€7r 

which equals the length of the first row in P, because those parts of a maximal 
path in n^jv which goes below the diagonal can be reflected in the diagonal to 
give a path in n^ 1 ^ without changing the sum ^ w (hj) since w(i,j) is symmetric. 
The same argument as above now gives 

F[F(N)<t] = -±j II Ihi-hjlflq^ 2 . 

N heN N l<i<j<N i=l 

max{hj}<t+N-l 

This corresponds to (3 = 1, 7 = 1 in theorem 2.2. It should be possible to analyze 
the asymptotics in this case analogously to GOE, see [TW2], to show that we can 
find constants a and b so that ¥[F(N) < aN + sbN 1 / 3 ] conerges to F±(t), the 
Tracy- Widom distribution for GOE. However it is not immediate to generalize the 
techniques of [TW2] , so this remains to be done. Note that again we can take the 
limit q — > 1 to get the case of exponentially distributed random variables. 

3. The Fredholm determinant 
From the identity (2.7) we have 

p N {t)=^ N (t + N-l), (3.1) 

where 

N 

^ N (s) = £ N [l[(l-Xs(h j ))}. (3.2) 

Here 

1 N 

Zm { n),n h ^ N 

K(N) = M(N) - N + 1, M(N) = [7JV] and Xs(t) is the indicator function for the 
interval (s, 00). We will take s in (3.2) to be an integer. 

Let M^' q (x), j = 0,1,... be the normalized orthogonal polynomials with 
respect to the weight w q K {x) on N, 

00 

£M^(s)MfV)<(*) = *y, 

x=0 
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and Mf' q {x) = Kjx j + ... with kj > 0. Set 



N-l 



K N (x,y) = M«>^x)M«>*(y)wJ c (x) 1 ' 2 w 9 K (y) 1 * 



3=0 



so that Kw(x,y) is a reproducing kernel on £ 2 (N). 

The polynomials M^' q are multiples of the standard Meixner polynomials, 
[NSU], [Ch], 

M«>«(x) = K —>-m^{x), 

where 

2 _ n!(n + K- 1)! 
n ~ (1-9)^^(^-1)!" 

The leading coefficient in m^ q is (^-) n and consequently 

_ 1 " 
d n \ q 

The Meixner polynomials have the generating function, [Ch], 

E m »^ = d--r(i-r" K (3.3) 

n=0 H 

The Christoffel-Darboux formula, [Sz], gives 

K N (x,y) = ^-i Mat^M^ -^^^-1^)^(^)1/2^(^1/2 

q m^W-ifa) - m N {y)m N - 1 {x) q , o/ 2 « , 0/2 ,« 4 n 
" (1 - g)^.! ^ k( J ^ ' ( ] 

where we have omitted the upper indices. Standard computations from random 
matrix theory, [Me], Ch. 5 and [TW2], show that ip^ can be written as a Fredholm 
determinant, 

N (-l) k 

^Pn(s) = k \ ^ det(iTjv(/ii,/ij))i<i,j<fc- (3-5) 

fc=0 ' he{s+l,s+2,...} fe 

The proof of theorem 1.2 is based on taking the appropriate limit in (3.5). 

The next lemma will allow us to compute the asymptotics of the right hand 
side of (3.5). 
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Lemma 3.1. Let b > be a constant and assume that pn — > oo as N — > oo. Sup- 
pose furthermore that Kn :NxB— > R, N >1, satisfies the following properties, 
(i) Let Mi > be a given constant. There is a constant C such that 

oo 

^2 K N (bN + p N r + m, bN + p N r + m) < C (3.6) 

m=l 

for all N > 1, r > -Mi. 
(j'i) Given e > 0, there is an L > so that 

oo 

^ K N (bN + p N L + m,bN + p N L + m) < e, (3.7) 

m=l 

for all N > 1. 

(j'ii) Let M > be a given constant. If A(£, rj) is the Airy kernel defined by (1.7), 
then 

lim p N K N (bN + p N £, bN + p N rj) = A(£, rj) (3.8) 

TV^oo 

uniformly for £, rj G [-M , M ]. 
(j'vj The matrix {Kjq{xi, Xj))^j =1 is positive definite for any xi, Xj £ [0, oo), k > 1 
Then, for each fixed t e R, 



lim V^f V det(K iV (6iV + p i vt + /i i ,6Ar + p iV t + ^))f . x = F(£), (3.9) 

where F(t) is given by (1.8). 
Proof: It follows from (iv) that 

k 

\det(K N (x i ,x j )) 1 <ij< k \ < Y\_K N (xj,Xj), (3.10) 

i=i 

see for example [HJ]. Consequently, 

(oo s k 

^2 K N (m,m) J . (3.11) 

where we have written a at = oiV + pArt. 
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Choose Mi so that \t\ < M 1 . Let e > be given. It follows from the estimates 
(3.6) and (3.11) that we can choose £ so that 



N 



( "i \ k /~^k 

E -IT- E det(K N (a N + h l7 a N + fe i ))f J . =1 | < ^ ¥ <e, (3.12) 



for all N > 1. Choose L so that (3.11) holds with L = L - M . Then, by the 
estimates (3.6), (3.7) and (3.10), 



- ^ det(KAr(aAT + ^i, a^v + ^))i<i,j< 

h€N fc he{[L i)PN ] c ) k J 
k 

< ^2 Y[K N (a N + hi,a N + hi) 

some hj>Loppf 
k k 

E i J /vyU'.y - /'M'.v • /', ) 

j = l /,, 6N fc i=l 



oo 



(OO \ fc— 1 s CO \ 

K N (a N + m,a N + m) \ I ^ K N {bN + Lp N + m,bN + Lp N + m) j 
m=l ' Vn=l ' 



< kC h ~ 1 e. 



(3.13) 



Denote the Fredholm determinant in the right hand side of (3.9) by D N (t). In- 
serting the estimates (3.12) and (3.13) into the formula (3.9) we obtain 



D »® ~ E { -TT- E det(/C*(<r + ^, a + ^L)) fc J_ 
fro fc! he[ £^ ^ ^ - - 



^(^^ + 1 )^ (1 + eC)e ' 



(3.14) 



where 



£aK£, ?/) = PNK N (bN + p N £, bN + p N rj). 



By assumption (iii), with Mq = Lq + Mi, we can chooose No so that if N > No, 
then 

\det()C N (a+ — ,a + — )) - det(A(a + — ,a + — ))| < 4r 
Pn Pn Pn Pn Lq 
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for all i,t/G [LqPn]. Thus, 



y (~i) fc v 

^ fe! ^ 

fc=o /ie[L pjv] fe 



< 



fc=0 



1_ / LqPat + 1 

fc! V l oPn 



det(lC N (t + ^±,t+^))- det(A(t +^,t+^)) 
Pn Pn Pn Pn 



e < C'e. 



Pn 



(3.15) 



Pn 



Combining the estimates (3.14) and (3.15) we find 

k=0 he[L oPN ] k 

The Airy kernel can be written, [TW1], 

/•OO 

A(x,y)= / Ai (x + s)Ai(y + s)ds 
Jo 

Using the formula, see for example [Ho], p. 214, 



< C"e. (3.16) 



(3.17) 



/ 1 f c 
Ai(x) = e-^ 32 -y_ 

valid for x > 0, we see that 



\Ai(x) \ < 



-e s 



2 r 3/2 

x , x > 0. 



2V7ra 1 /4 

This estimate can be used to show that the Airy kernel satisfies (i) and (ii) above. 
Since the matrix (A(£i, £j))i<i,j<k is positive definite, we can use the same argu- 
ment as above to show that 



£/ -E 



k=Q J[t,oo)>* k=0 J[t,L ]k 



^det(A(^))* J=1 ^ 



< e 



(3.18) 



provided £ and L are sufficiently large. From (3.17) we see that choosing Ni > N 
large enough we have 



det(A(^,^-))i<i,j<fc^ 



< C""e 



(3.19) 



for all N > Nx. If we combine the estimates (3.18) and (3.19) we have proved the 
lemma. 

To apply this lemma to the Meixner kernel (3.4) we need 
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Lemma 3.2. The Meixner kernel satisfies the properties (i) to (iv) in lemma 3.1 
with b = 6(7, q) as before and pn = crN 1 ^ 3 , where a is given by (1.11). 

This lemma will be proved in section 5. We can now combine (3.1), (3.5) and 
(3.9) to get 

lim p N ((b - 1)N + aN 1/3 t) = F(t), (3.20) 
which is (1.10) and theorem 1.2 is proved. 

4. Proof of the large deviation theorem 

In this section we will prove theorem 2.2. Set 

K N y= ^ k V {Xi,Xj). 

By adding a constant C to Vn, which does not alter the problem we can, by 
assumption (ii) on Vn, assume that 

V N (t) - log(t 2 + 1) > £log(t 2 + 1) (4.1) 

for all t > 0. Since \t - s\ 2 < (t 2 + l)(s 2 + 1), this implies 

M-l 

-K M ,v N {x) < -£(M - 1) Y, lo g(! + *?) ( 4 - 2 ) 

i=i 

for all x £ [0, 00) M . Note that 

N-l N-l 

Y log|xj - Xk\ - N Y v n(xj) = -K N -i,v N (x) - Y v n(xj). (4.3) 

l<j^k<N-l j = l j = l 

The next lemma is analogous to lemma 4.2 in [Jo]. 

Lemma 4.1. Let {sn} be a sequence in [0, 00) such that sn — > s > as N — > 00, 
or sn = 00. Set, for a given a > 0, 

fWOO = {x E A JV ( S ) 7V - 1 ; l^K N -i,v N {x) < Fy + a}. 

Let < A < 1 and let on G A at, N > 1, be a sequence converging to a > 0. Define 
a probability measure on Atv^tv)^ -1 by 

P n '-In(^sn) = , ON E II I^-^I^Qjv-i^Cx), (4.4) 

z n-i,nK s n) x€ q j=i 
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where Z^°_i n( Sn ) is a normalization constant. (E^li N [-', sn] denotes the corre- 
sponding expectation and if sn = oo or A = o~n = we omit them in the notation.) 
Fix n > 0. Then there is an Ni such that for all a > and N > Ni, 



P^Zl N (n N ,r,+a(s N ) c ;s N ) < e~-^ N \ (4.5) 
PROOF: We first prove the following claim. 

Claim 4.2. Let G A^v, on — > a as N — > oo and s G (0, oo]. For each N > 2 
we can choose (xf , . . . , £jv_i) G Aa^s)^ -1 so that 

JV-l iV-l 

^ E lQ g i-f -^r 1 + ^E ^v«) - ^ E lQ g ^ - < i - ^ 

l<j/fc<Af-l i=l j=l 

(4.6) 

as A" — > oo. 

To see this set 

7 r j/N k 

y% = max{^ ; j G N and y <^(t)eft < -}. 

If ^ o-jv for = 1, . . . , N-l, we set x% = y% . If y^ Q = cr N , we set x% = y£ for 
k < k and x% = y% + 1/N for k = k , . . . , N - 1. Using the fact that < (j) 8 v < 1 
it is not difficult to see that < < . . . < x^_ 1 < L for all A" and some fixed 
L. Furthermore 



I N ~ 1 

/V^T E ~* <t>v( x ) dx ( 4 - 7 ) 
k=l 



weakly as A" — > oo. The property (iii) in the assumptions on Vjv implies 



.7 = 1 70 



Clearly, 

^ E lQ g - *f r 1 < ^ e lo § y = lo s (^3i)i ' ( 4 - 9 ) 



which — > as A" — > oo. Also, since o/v — > a and the belong to a bounded set, 
we get a bound in the other direction which goes to as A" — > oo. Given M > 1, 
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set f M (t) = min{log|t| -1 ,logM}. Write 



N 2 
+ 



1 

N 2 



Qog\x?-xS\-l-fM{x?-xS)). 



(4.10) 



3+k 
-*>k\<l/M 



The absolute value of the second sum in the right hand side of (4.10) is 



< 



N 



N 2 



l<\j-k\<N/M 
\j\,\k\<LN 



,logM 
M ' 



Thus, using the weak convergence (4.7) and then letting M — > oo we obtain 

2 poo poo 

i im E^sl^f -^l" 1 = / / l °s\ x -y\~ 1( t>v( x )<f>v(y) dxd y> 
N ^°° N ttk Jo Jo 

which together with (4.8) and (4.9) proves the claim. 

We turn now to the proof of lemma 4.1. Let e > be given. We want to 
estimate Z^°^ N from below. Choose No so that s^v > s — e ii N > Nq. Then 

z n-i,n( s n) > z n-i,n( s - e )> 
if iV > N . Choose (x% )^ = " 1 1 C A N (s - e) as in the claim. Clearly, 



jlp log Z%Zl )N (s N ) >-- 



\ x j x fc I 



j^k 



N-l 



N-l 



JVi 



and consequently, by Claim 4.2, 



iMf ^ lOgZ^^) > ~^Fy~ 



Since F* _£ \ F° as e -> 0+, 



Ij^mf j^logZ^_l N (s N ) > -'-F^. 



(4.11) 
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Thus, given 5 > 0, we can choose N(S) so that if A" > N(S) ,then 

j^logZ x ^_l N (s N ) > -£{Ft + S). (4.12) 
It follows from (4.2) with M = N - 1 and (4.3), that for any < p < 1/2, 

-^-f,Af(^A r ,'7+«(^A r ) C ; s n) 

N-l 

r ~i N 

< e ^(F^+5)-§(l-p)(F^+ v +a)N 2 (t 2 + l)-^^-^ (1 + a%) X0/2 

JEA N 

if A" is sufficiently large (independent of a > 0). Note that 5 + pFy — rj < if we 
choose 5 = rj/2 and p sufficiently small. This completes the proof. 

This lemma can be used to prove 
Corollary 4.3. For any s e (1, oo], 

Furthermore F^ - F v > if s < b v . 

PROOF: The lower limit follows by taking A = cr^ = in (4.11) (replacing N — 1 
by N does not modify the argument above in any essential way). Given < p < 1, 
we can use (4.2) with M = N and the continuity of expKN,v N to see that 

Z N (s)= Yl e~^ KN ^ {x) ~^ N ^ VN{xj) 
xeA N (s) N 

< sup e -^- P )K N ,v N {x) J- e-l^-^E,- 10 ^*?) 

< e -|(i-p)K N ,v iV (y JV )+C7V j ^ 4 _ 14 ^ 

if A" is sufficiently large, where = (yi , ■ ■ ■ ,Vn) e &n{s) n ■ Clearly, yj^ 7^ y^ 
if j ^ fc. Set X N = N' 1 J2j V- Zt follows from ( 4 - 12 ), with A = a = and N - 1 
replaced by A 7 ", that N~ 2 log Z N (s) > -(3{F^ + S)/2 for N > N(5), so (4.2) and 
(4.14) yield 

\og{\ +t 2 )d\ N {t) < C. 
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Thus {Aat}^ =1 is tight. Pick a subsequence that gives the upper limit of 
N~ 2 log Zn(s), and a further subsequence so that Ajv^ converges weakly to v = 
ifjdx. The measyre v has to be absolutely continuous with density satisfying 
< ip < 1 because of the definition of A at. Using (4.1) and \t—s\ < \ft 2 + 1a/s 2 + 1 
we see that ky N (t, s) > 0. Set, for given M > 0, ky N (t, s) = mm(ky N (t, s), M) and 
choose 4> T {t) continuous so that < T < 1, T (t) = 1 if \t\ < T, = if \t\ > T+ 1 
and (j) T (t) < (p T '(t) if T < V . Then, k VN (t,s) > (h(t)(h(s)k^ N (t, s) and using 
the estimate (4.14) we get 

— J log ZAr.(s) 

and thus, letting j — > oo, M — > oo, T — > oo and p — > 0+ in that order, we obtain 

-^F v ^limmf-^logZ^s) ^limsup-^logZjvCs) < -^vM- 

Thus [-0] < -Py and ip £ A s , so we must have ip = 4> v . 

Assume that F v < F v and s < by. Then ly^y] < ^v[0v] an d consequently 
(p v = <py by the uniqueness of the minimizing measure. This contradicts the 
definition of by. The corollary is proved. 

Note that by (2.13) corollary 4.3 implies (2.16) so we have proved the first 
part of theorem 2.2. Before turning to the proof of the second part we need one 
more consequence of lemma 4.1. 

Corollary 4.4. Let {sn} be as in lemma 4.1 and assume that f : [0, cr + e] — > R, 
e > 0, is continuous, or f : [0, oo) — > K is continuous and bounded in case sn = oo. 
Then ^ 

lim 1 ]QgE%™ N [e^U s N ] = [°° f{t)Mt)dt. (4.15) 

Furthermore let 

i=l 

(8 t , s is Kronecker's delta), be the 1-dimensional marginal distribution of the prob- 
ability measure (4.4) (with sn = oo). Then for each < y < 1: 
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(ii) if S t is the Dirac measure at t, then X^teAjv w jv-T jv(*)^* converges weakly to 
4>v{t)dt as N — > oo. 

PROOF: We can prove (4.15) using lemma 4.1 in exactly the same way as in the 
proof of (2.5) on p. 194 in [Jo], see also [De]. The weak limit (ii) is a direct 
consequence of (4.15), see [De]. Note that the limit does not depend on y since 
the factor \ a N ~ x i\ vl3 does not affect the leading asymptotics. 

In the expectation (4.16) all the Xi :s have to be different, otherwise the prob- 
ability is zero, and consequently the expectation is < 1, which proves (i). That 

(iii) holds follows from the presence of the factor \ a N — Xi\ vl3 - The corollary 
is proved. 

We turn now to the proof of the upper-tail limit. Note that 

M-l 

Qm,n(x) = e -^^(xM) "Q \ Xm _ Xi fQ M _ hN ( x '), (4.17) 

i=i 

where x' = (xi, . . . , xm-i)- Using this identity we see that 

(t) = Ml Qm,n(x) 

Xl<...<XM<t 

M-l 

seAjv(t) x6A JV (s) M - 1 i=l 

If we define 

H M -1,n(s) = p-r 1 f \S - XifQ M -l,N(x) 

Zm"_i n(s) — 

i,JV V ' x6A J v(s) M " 1 i=l 

this can be written 

Z M , N (t) = M e-^ VN ^Z M -i,N(s)H M -i,N(s), (4.18) 

s6Ajv(t) 



or 



F M ,N(t) = MZ y M ^ N e-^ v ^F M ^ N (s)H M ^ N (s). (4.19) 

Zm > n settit) 
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This is the main formula to be used in the proof of (2.17). We will need two choices 
of M, namely M = N and M = N — 1. They are handled completely analogously 
and we will consider only the case M = N. 

Write Ajv(t, s) = An n (t, s) for any < t < s < oo and Aat(£)* = Ajv(t, oo). 
If we let t — > oo in (4.19) and then subtract (4.19) from the limiting equality, we 
get 

l-F N (t) = N i N - 1 > N J2 e-^ VNi8) FN-iAs)H N -iAs)- (4-20) 

Z/ AT M 

s€A N (t)* 

Set 

1 f°° 

$v = F v -- J V{s)(j) V {s)ds, 
From the variational relations for <py (t) it follows that 

/ log|&y-s| - Vv(s)ds + ^(fty) = $v- (4-21) 



Lemma 4.5. We have 



limsup-Log ^- 1 '" <P$ V . (4.22) 



Proof: By (4.17) we have 

^ e -*v»(., B [JJ |, 
^n-i,n rr 

AT-l 

>e-^ v ^)E N _ hN [U \r-Xi\f] (4.23) 

i=i 

for any r E A N . One difficulty in estimating the right hand side in (4.23) comes 
from the fact that, due to the discrete nature of the problem the integrand could, 
apriori, be zero for many y :s with high probability. Note that we define y = 
for any y > 0. Let ip a (t) = 1 if t ^ s and ip s (s) = 0. 
Consider 

1 N-l 

S) = — logE N _ 1)N [ JJ |s - Xi| tf ^a(^i)]- 
i=l 
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Then, 

1 N-l 
f N (0+;s)= lim f N (y;s) = —\ogE N - 1)N [T[ip a (xi)] 

i=l 

= — log Pat_ i, at [all x % ^ s]. (4.24) 
Let e > be given and write Bjv(e) = An (by + e, by + 2e). Now, 
^2 Pv-i,iv[all x i ^ s] > P N -i,n[ (J {all Xi ^ s}] 

s6Bjv(e) seljv(e) 

= 1-P N _ 1 , N [ P| {one Xi = s}]. (4.25) 

s6Bjv(e) 

Take g : [0, oo) — > [0, oo) continuous such that g(,s) = 1 if by + e < s < by + 2e 
and g(s) = if < s < by or s > by + 3e. Then, 

e^Pjv-i,*! f| {onex^^^iv-!^^^]^^ 2 (4.26) 

s6Bjv(e) 

for all sufficiently large N. The first inequality follows from the definitions whereas 
the second follows from corollary 4.4, (4.15). Combining (4.25) and (4.26) we see 
that 

max PAr_ ljA r[all x % ^ s] > — (4.27) 

s6B JV (e) 2N 

for all sufficiently large N. Hence, by (4.24) and (4.27) we can choose = 
cr7v(e) G Iw(e) so that 

lim /jv(0+;<7Ar) = 0. (4.28) 

Take r = a N in (4.23). Then 

1 t Z NN 13 . . 

— log > -- VW(^jv) + /iv(l; ^jv) 

iv Zn-i,n * 

= ~V N (<r N ) + f N (0+;<T N )+pJ f' N (y;a N )dy. (4.29) 

We can pick a subsequence {Nj} which gives liminfAr^oo log Z ^ N '^ and such 
that a Nj (e) -> a(e) G [6 V + e, 6 V + 2e]. Then, by (4.28) and (4.29), 

liminfllog^^- > -^(a(e))+/31iminf / f' Nj (y, a Nj )dy. (4.30) 
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Now, 



N-l 



f N (v> °n) = En-iMjJ log \ aN ~ Xi 



N-l 



N 



Hence, by corollary 4.4 (i) and (iii), 



\og\a N -t\u^_l N (t). 



teA N 



N 



>2-^X>g^> -2 

i=l 



and consequently, by Fatou's lemma, 



liminf / f' N (y; a N )dy > / liminf f' N (y; a N )dy. (4.31) 
Jo Jo 



Given 5 > 0, small, and M > set 



logM, if |tj > Af 
fM,s(t)=< l °s\tl H5<\t\<M 
log 5, if \t\<S. 



By corollary 4.4 (i) and (iii) we have 



(min(logM, log \a N - t\) - f M , s (<r N - t))u y ^_ N 1N (t) 

teA N 



* E 

teA N ■ 0<\t-a N \<8 

2N 



log 



iV-1 - iV-1 ^ 5 k 

k=i 



< 



-S. 



N-l 

Also, if | a n — o~e | < 5, which is true if iV is large enough 



|/m, 5 (K -*|) -f M , s (We) -t\)\ < Slog-. 

Since log \(Jn — t\ > min(logM, log \ctn — t\) and M,S are arbitrary it follows from 
corollary 4.4, (ii) that 

POO 

liminf f' N (y; a N ) > / log \a(e) - t\(j) V (t)dt. 

3^c° 3 Jo 
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Together with (4.30) and (4.31) this gives 

limmf-log-^^>-^(a(e))+/3 / log \*(e) - t\$ v (t)dt. 
We can pick a sequence ej — > such that <r(ej) — > by and using (4.24) we obtain 

i- • r 1 ! ^ N,N 

hmmf — log — - — > 

N^oo iV Zjv— 1,JV 

and the lemma is proved. 

Given 5 > we can use lemma 4.5 to find N (S) so that 

Zn,n 

if N > N (6). Since Fiv-i,Ar(s) < 1 we can combine (4.20) and (4.32) to get the 
estimate 

1 - F^t) < iVe^*^ ]T e-^^W^_ liJV (s). (4.33) 

s6Ajv(t)* 

We have 

N-l 



H N - ljN (s) = E N _ ljN [Y[ \s - Xif; s] 

i=l 

< (1 + S 2 )^ (Ar - 1) ^°-i,Ar[ ^(l + ^)^ /2 ; s] < e CN (l + sY N '\ 



i=i 

where the last inequality is proved, using lemma 4.1, just as (4.25) in [Jo]. Together 
with (4.1) this gives 

e-¥v»i>)H N - ltN {8) < e CN-^io g( iW) : (4 _ 34) 
Hence, given a constant D > 0, there is a constant d > such that 

e N Pi * v+s) J2 e- N ^^/ 2 H N _ ljN (s)<e- ND . (4.35) 

seA N (d)* 

For t > s we define 

H N -i >N (t, s) = — y~] jj \t-XifQ N -i,N(x). 
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Clearly, 

H N -i jN (s) = H N - ljN (s, s) < H N - ljN (t, s) (4.36) 
if t > s. Combining the estimates (4.33), (4.35) and (4.36) we obtain 

l-F N (t)<Ne- ND +Ne N ^+V £ e~^ v ^ H N _ hN (s + e, s) (4.37) 

for any e > 0. Let sat G Ajv(t, d) be the s which gives the largest term in the sum 
in (4.37). Then 

1 - F N (t) < Ne~ ND + N 2 (d- t)e N ^ v+s -^ VN ^H N _ ljN (s N + e, s N ). (4.38) 

Choose a sequence which gives the upper limit of iV _1 log(l — i 7 V(t)) and such 
that — > a G [£, d]. We would like to prove that 

lim — logH N ljN (s Nj + e, Sjv .) = -/3 / log |a + e - t\(j)^(t)dt. (4.39) 
./ • * iVj ./ 

We will write N instead of Nj for simplicity. Looking at the definition of 
HN-i,N(t, s), we see that we are interested in the limit of 

1 i 7-i r * log |sjv+e— Xi\ n 

as iV — > oo, — > a. Since 

| log \s N + e — xA - log \a + e — xAl = | log |H — — — — 1| < C-^ — — , (4.40) 

a + e — a + e — 

where C is a numerical constant, and sn < c + e/2 for iV large enough, the limit 
(4.39) follows from corollary 4.4. 

If £ > by, then 0^ = 0^, since cr > t, and combining (4.38) and (4.39) yields 

limsup-J-log(l - F N (t)) 

< max{-.D, / 0$y + <$- ^V(a) - 13 J log \a + e - t\~ l (j) V (t)dt}. (4.41) 

Note that a could depend on e and d. Pick a sequence e = e-,- — > 0+ and then a 
subsequence so that cr(e Jfe ) — > r G [t, d]. Then, since D and 5 are arbitrary, we get 

limsup — log(l-FAr(£))</?($ y -inf [ k v (r, s)(f) V (s)ds) (4.42) 
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and we have proved one half of (2.17). 

We now turn to the lower limit. If we start with M = N — 1 instead of N then 
(4.42) holds with F N _i replaced by F N _i iN (t). By assumption the right hand side 
of (4.42) is negative for all t > by. Hence, if t > by, we see that 

F N -!, N (t) > 1/2 (4.43) 

for all sufficiently large N. Note that, if t > s, then 

H N - hN (t) > ZN - hN [ S } H N ^ N (t, s) > F N - hN (s)H N _ hN (t, s). (4.44) 

^N-l,N{t) 

The function f(r) = f ky(r, s)(f)y(s)ds is continuous on [t, oo) and f(r) — > oo 
as r — > oo, so it assumes its minimum in [t, oo) at some point tq > t. Let e > 0. 
Pick sn £ A;v(t)* such that sn \ To + e. Then, picking one term in the sum 

e-^ v ^F N _ hN (s)H N _ hN (s) 

seA N (t)* 

> e -^ v »( s ^F N _ ljN (T ) 2 H N _ hN (s N , s N - e). 



If we use the limit (4.39), the estimate (4.43) with s = To, and let e — > 0+, we see 
that 



liminfi-log V e-^^^F JV _ 1)J v( S )^_ 1)JV ( S ) 

N^oo iV — * 

s6Ajv(t)* 

> -§^(to) - log |t - tl-^y^dt. (4.45) 



To complete the proof we need 

LEMMA 4.6. For any Vn satisfying the conditions (i) - (in), 



lim inf — log Zn ~^ n > p$ v . (4.46) 



Proof: If we let t — > oo in (4.19), we see that, e > 0, 



Zn-i,n 



< N J2 e- — VN ^F N _ 1 Ms)H N -i,N(s) 

seA N (b v —e) 

+ N e-^ v »^H N _ hN (s), (4.47) 



s€A N (b v -e)* 
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since F N _i ;N (s) < 1. By adjusting the constant C we see that (4.34) holds for all 
s G Ajy, so the first sum in the right hand side of (4.47) is 



< e CN F N ^ N {b v - e) e-^^d+s 2 ) < e CN-§L( bv - £ )N- 



s6Ajv 



for all sufficiently large N by the first part of theorem 2.2. (Replacing F^(t) by 
FN-i,N(t) does not make any difference.) Since L(by — e) > if e > 0, the first 
part of the right hand side of (4.47) is negligible. 

The same argument that lead us from (4.33) to (4.42) allows us to treat the 
second term in the right hand side of (4.47) and obtain 

r 1 , Z N ,N 

iim sup — log 



< max{-L>, ~V{a) - (3 J log \a + r] - t\~ l (j)^{t)dt}. (4.48) 

where a G [by — e, d] , r\ > 2e, D > are given. Take e = ej — > 0+ so that 
cr(ej) — > t G [by, d] . Note that 4>p ej \t)dt converges weakly to <fiy(t)dt = 4> v {t)dt. 
Using an inequality like (4.40) we get 

,. 1 , z N ^ N 

hm sup — log 



N^oo N Zjv-1,JV 

< max{-L>, ~V(t) - (3 J log \t + rj - t\- l (j) V (t)dt}. (4.49) 
We can now repeat the argument that lead from (4.41) to (4.42) and obtain 

1 N (3 f f 

limsup — log — < - / V(s)(j) V (s)ds - (3 inf / k v (r, s)(j) V (s)ds < -(3§ v , 

N^oo Zjv-1,JV * J T > b v J 

since J kv(r, s)(pv(s)ds > F v if r > fry. The lemma is proved. 
Combining (4.20), (4.45) and lemma 4.6, we see that 

liming log(l--M*)) 

AT— >oo iV 

> (3(F V - J k v (r ,s)(f) V (s)ds) = (3(F V - inf J k v (r, s)(f) V (s)ds), 



by the choice of tq. This completes the proof of theorem 2.2. 
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5. ASYMPTOTICS FOR THE MEIXNER KERNEL 

This section is devoted to the proof of lemma 3.2, which is based on establish- 
ing the appropriate asymptotics of the Meixner polynomials. See [Go] and [JW] 
for some results on the asymptotics of Meixner polynomials. 

From (3.3) we obtain, x e R, 



K ~n+l 



sin the {y/l) n+K n\ f r 



y/7-t/y/q 



dt 

(5.1) 



~ y/&) K t n+1 



where T r is the circle \z\ = r, < r < yj^/q; if < r < ^fq the second integral 
should be omitted. Let b = (1 + v /tq) 2 /(1 — q) as before, let a be given by (1.11) 
and set 

(v/7 + v^) 2 



a = b + 7 — 1 



Set 



and 



t(z) = 
s(z) = 

A N (x) = 



1-q 

iq + z\ ( y/q 



y/l + ^/qz 
b x (x + K- l)\M~f K+N fq 



x x + K x\(N + K - 2)! 1-q V 7' 
For < r < a/tA? we define 

^n^; #) = if r > v /t1/ and 

^fo</) = (-l) n+ * +1 / |t(-r)r a (-r)^(-r)— pj. (5.4) 

The powers are defined by taking the prinipal branch of the logarithm. 

The Meixner kernel (3.4) can now be written, for x, y integers (which is the 
case we need), 

K N (x, y) = V A N (x)A N (y) 9l) ° N ^ ^ = 9l) (5.5) 

x-y 
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if x 7^ ?/, and 

K N (x,x) = A A r(x)[ J D 7 v(a: - 1; g 3 )D N (x; g 2 ) - D N (x; g 1 )D N (x - 1; g 4 ) 

+ F N (x;g 1 )D N (x;g 2 ) - F N (x; g 2 )D N (x; gj], (5.6) 

where gt(z) = 1, g 2 {z) = z - 1, g 3 (z) = t(z) \ogt(z) and g A (z) = g 2 (z)g 3 (z). The 
functions gi(z) are bounded for \z\ < 1. 

Write x = Nb + y and K = [7JV] - N + 1 = iV^ - 1) = #(7 - 1) + oo N , 
<u> N <1. 

Lemma 5.1. If x = Nb + ^criV 1 / 3 and M > is a given constant, there are 
constants 01(5,7) an d c 2 (q,"f), such that 

^A N (x)<c 1 (q^)e c ^^ N ~ V3 (5.7) 
for all £ > —Mo. Furthermore, 

lim ±=A N (x) = 7 ^ (5.8) 

uniformiy for |£| < M . 

PROOF: By Stirling's formula 

(x + K) x+K N N b x g+jV (iV + K)(iV + K-l) 
N{X) ~ x*{N + K) N + K a*+ Kl x + K 



Write a,N = b + 7^ — 1. Then, 



(x + K) x + K N N F ^k+n = (Nb\ x (x + K\ x+K (a N \ x+K ( 7 ^ V+A 



x x (N + K) N + K a x+K 1 V x J V / \o/ \7n 

(5.10) 

If we write u — Na^ and t> = Nb < u. Then 

/ a 7" 7 \ x / , T r\ x+K / \ u+y / \ — v — y 

fNb\ fx + K\ , V\ , V\ ■ _g(y) 



Since g(0) = g'(0) = and g"(t) = (v - u)(u + t)' 1 ^ + t)' 1 < 0, we have 
exp g(t) < 1 if £ > 0. If -M < f < M , then 



|</(t)| = | / (t - < CN~ 

Jo 



1/3 
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Furthermore 

\ x+K 

a N \ = u N+ o(dN- 2 / 3 )+o(l) 



a 
and 

K+N 



_L] =e -co N +o(l)_ 



.In. 

Inserting these estimates into (5.10) we obtain 



(x + K) x+K N N b x ^ K+N ^ Ce c iN -^ 



x x (N + K) N + K a x + K 
for £ > -M and 

(x + Ky^N^ K+N 
iv" oo x x (N + K) N+K a x+K 

uniformly for |£| < Mq. By (5.9) this proves (5.7) and (5.8). The lemma is proved. 

Set 

u(z) = b log( v /7Q + z) - a log(y/7 + Jqz) - log z 

so that 

D r N (x;g) = — [ e N(u( z )-u(i))+yio g t( z )+co N io g s( z ) g ^d^_ (5 _ n) 
2tt J r iz 



Now, 



u'(z) = -p(l-z) 



2 



(1 z \3 Vv z2 + (y/Q + Vl + Qy/^ z + VQ + Vl + Qy/l + lV<i 

^(^ + v / 7?)(v / 7 + ^) 



where 

p = 



Hence we can write 

1 



u{z) - u(l) = |p(l - zf + p(l - zfv{z), (5.12) 



where one verifies that \v{z)\ < 28/27 if \z — 1| < 1/4. 
By taking absolute values in (5.3) we obtain 

X I 2 „K(^ „\K pit 
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where 

K , o s 

/( T ) = 2 lo S(79 + r +2 v /7Qrr)H - — log(7 + qr +2 v ^rr). 

Write r = 1 — S, < 5 < 1. A computation shows that /'(r) > if (say) 

„ > V + g + VT^ (6 . 14) 

which covers all the y :s we are interested in. Thus, if (5.14) is fullfilled, then 

\D r N (x;g)\ < CeMN(u(l - d) - u(l)) + ylog*(l - d)). (5.15) 



By (5.12), 

if < 5 < 1/4. Now, 



u(l-S)-u(l)< P 5 s ( 1 -5^)<jpS 3 (5.16) 



log*(l -5)= log(l (1 J )v ^ -6) 

<-p(l-g)-L<S, 
v 7 ^ 



and consequently it follows from (5.15) and (5.16) that, if y > 0, then 

2 /V 1 
1^(0^)1 < C[exp[— p5 3 - p(l - «)— <5y]. (5-17) 

V79 



Recall that y = aN 1 / 3 ^ with a given by (1.11). Note that a = (1 - q)~ x ^qp~ 2/i ■ 
Choose 5 = (piV)" 1 / 3 ^ if £ < (A^p) 2 / 3 /16 and 5 = 1/4, if f > (iVp) 2 / 3 /16. 
Inserting this into (5.17) gives 

\D r N (x;g)\ ^Cexpf-imi^^^iVp) 1 ^], (5 . 18) 

for £ > 0. 

Let e G [0, it] and set 

1 f , .-()> . / .-()>«. / .-a, ^ 



I>1 = 1 LJ 9(re ie )t(re w Ys(re^) K 
I'{ = D r N (x-g)-I[ 
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By the same argument that was used for (5.13) above, we see that if y satisfies 
(5.14), then 



< C\t(re ie )\ x \s(re 



ie\\K 



,N 



< C exp[NRe (u(re ie ) - u(l)) + y log \t(re l 



(5.19) 



Next, we consider FJ v (x;g), ^Pfq < r < 1. Taking absolute values in (5.4) 
yields 



\F r N (x;g)\<C 



<yq-T 



79 + 1 



x+K 



dr 

-N+l 



(5.20) 



The integrand in (5.20) is a increasing function of r for all x that we are con- 
sidering. The monotonicity argument used for (5.13) now shows that, if (5.14) is 
fulfilled, then 

\Fr N ( X ;g)\<C\t(-rT\s(-r)\ K ± 



K 



-N 



< C\t{re ie )\ x \s{re ie )\ 

< Cexp[NRe(u(re* e ) - u(l)) + ylog\t(re ie )\], (5.21) 

where the last inequality is the same as in (5.19). If we take e = 0, we get the 
same right hand side as in (5.15) and hence we obtain the same estimates, i. e. 

I F r N (x; g) | < C exp [ - i min( yft, I (ATp) . 



Combining this with (5.6), (5.7) and (5.18) yields 

\K N (x,x)\<CNexp[-] min( ^ - (Np) Wft] 



(5.22) 



for any £ > 0; x an integer. 

Consider now f e [-M , (pN) 1 / 6 ]. Take e = (pN)' 1 / 4 , 5 = ^(pN)- 1 / 3 < 
(piV) -1 / 4 , where r\ > will be chosen below. By (5.12), we have 



J9\3 



I[ = ^ g((l - 5y d ) expjiV [ip(l - (1 - 6)e* 

+ p(l - (1 - dy e ) 4 v((l - 5)e ie )\ + ylogf((l - 5)e ie ) 

+ u N \ogs((l-5)e ie )}d6. (5.23) 
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We make the change of variables 9 = u(pN)- 1 / 3 . For < t] < (pN) 1 / 12 , \9\ < e, 
we have 

l - P {i - (i - 5y e f + P (i - (i - 5yy v «i - 5y e ) = ±( v - tuf + r u (5.24) 

where R\ — > uniformly as N — > oo. Furthermore, if £ G [—Mo, (piV) 1 / 6 ], then 

y log*((l - 5)e ld ) = (- V + iu)i + R 2 , (5.25) 

where -R2 — ► uniformly as iV — > oo. 

Suppose ^(1) = 0, j = 0, 1 but gW(l) ^ 0, so that 

g((l - 5Y e ) = IgM (l)p- £ / 3 (-V + ™Y + ■■■■ (5-26) 
We now have all the estimates we need. Let rj = y 7 ^ if £ > Mq and rj = 1 if 

ici<m . 

By (5.12) and (5.24) we obtain 

Re Nu((l - 5y d ) = \rf - r]iv 2 + R x 
and hence, if £ G [-M , (piV) 1/6 ], e = u(pN)- 1 / 3 with u = (pA^) 1 / 12 ,(5.19) yields, 

< Cexp[^ 3 - V (PN) 1/6 -vt + Rs] 
<]^IW-p[-^|3/ 2 ]. (5.27) 

Similarly, by (5.21), for £ G [-M , (pN) 1 / 6 }, 

W\ < j^j-^M-lKn- (5.29) 
The dominated convergence theorem gives 

n -(£+l)/S i roo ■ 

Km N( £+1 V 3 I[ = t— d^H 1 )^ J (-V + ^yexpl-iu + irif + i^tv + i^duj 

-(e+i)/3 

= P — X ^(l)AiW(O, (5.30) 

uniformly for \i\ < M . Observe that g x {\) = 1,# 2 (1) = but g' 2 (l) = 1, g 3 (l) = 
but g' 3 (l) = p(l-q)( iq )- 1 ' 2 and g 4 (l) = g' 4 (l) = but gjf(l) = 2p{l-q){ 1 q)- 1 ' 2 . 
Combining (5.27) and (5.29) we obtain 

\D r N (x;g)\ < ^TTW ex P[-^! 3/2 ]' ( 5 - 31 ) 
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for £ G [-M , (pN) 1 / 6 }. The estimate (5.27) and the limit (5.30) give 

lim N^D r N {x- 9l ) = p-^Ai (£), (5.32a) 

N^oo 

lim N 2 / 3 D r N {x;g 2 ) = p- 2/3 Ai'(£), (5.32b) 



Jim N*'*D* N {x;g 3 ) = ^L^Ai'fc), (5.32c) 



lim ND r N (x; 94 ) = ^— ^Ai"(£), (5.32d) 

iv— ><x> ylQ 



and 



We can now use (5.22), (5.28), (5.31) and (5.32) in (5.5) and (5.6) to prove (3.6), 
(3.7) and (3.8) for the Meixner kernel. The lemma is proved. 

6. The equilibrium measure 

The equilibrium measure (f>y(t)dt satisfies certain variational conditions. 
Proposition 6.1. Assume that 4> e A s satisfies 

(i) ^k V (t,T)(f>(T)dT>Xif(f>(t) = 0, 

(ii) /; k v (t, T)<t>(r)dT < A if (P(t) = 1, 
(Hi) f* k v {t, t)0(t)oV = A if < 0(£) < 1, 
for some A (which = F v ). Then = 0y . 

We will not prove this here, see [LL] for a very similar result. The way to 
compute 4>v is to seek a candidate solution and then verify that satisfies the 
variational conditions. In a region where < 0(£) < 1 we can differentiate (iii) 
and obtain 

1- 



dr = --V'(t). (6.1) 
o r ~t 2 

Since V 1,q is convex the support of 0y is a single interval. If we consider the 
variational problem without the constraint < 1, and this problem has a solution 
ipo such that < ip < 1, then this ip is the solution we are seeking. This is the 
case when 7 > 1/q, and then [ay, Oy] = [a, 6] and 

'r/v a<t<b. (6.2) 



r - t 2 



We must have 0(6) = and 0(a) bounded (0(a) = if 7 > 1/g). 
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If the solution ipo(t) > 1 in some interval, e.g. tpo(t) > 1 in [0, a ) but 
< ipo(t) < 1 in (oo, bo), we make an ansatz that <f>(t) = 1 in [0, a] and < <f>(t) < 1 
in (a, b) for some a, b, [ay, by] = [0,6]. This is the situation when 7 < 1/q. By 
(6-1), 

and 4>(a) = 1, 4>{b) = 0. By making the substitution x = 2(t — a)/c — 1, y = 
2(t — a)/c — 1, c = b — a, in (6.2) and (6.3) we get an equation of the form 

1 ' ' U) -dx = f(y), -1<x<1, (6.4) 



K J-l 



x-y 

with some /. This equation has the general solution, [Tr], 



nVl-x 2 J-i y-x nVl-x 2 

where C is an arbitrary constant. In this way we obtain (2.19) and (2.20). 

The equation (2.21) is obtained by substituting (2.19) or (2.20) into (2.15) 
(the infimum is assumed for r = t). Consider the case 7 > 1/q, the other case is 
similar. Then, with t = a + c(x + l)/2, 



J(t) 



J J\s)ds=^J J'(a + c(y+l)/2)dy 



and 



g(y) = J'(a + c(y + l)/2) = °- f ^dx + \v'{a + c(y + l)/2) 

^ J -i x — y z 

f 1 11 

/ log \y - x\v'(x)dx + -[log - - log(y + B) + log(y + D)]. 
J-i 1 Q 



c 
2 



Now, 



and 



where 



VD 2 - 1 VB 2 -1 
x+D x+B 



f 1 11 

J log \y - x\v'(x)dx = -F(y, D) - -F(y, B), 

1 f 1 V 1 R 2 — 1 
F(y, R) = - - — , log \y - x\dx. 

7T J_! [ X + R)Vl - X 2 
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Note that 

d „, VR 2 - 1 r 1 1 

-F(y,R) = ■ - [-=— + 



Using these formulas we see that gr(— 1) = and hence 

c f x c f x 

J ( t ) = ^J 9(y)dy=-j (x-y)g'(y)dy 

[x-y){- 



which gives (2.21). 

If /(y) = (7 - q)(v + B)- 1 + (1 - g 7 )(y + -D) -1 , then f(y) > for all y > 1 
and a = infi< tf <i/ c /(y) > 0. Thus for < 5 < 1, by (2.21), 

J(b + 5)>^ f 1+5/C (l---y) . ^ > Cl ^ 2 , 
8 v /<ry y x c + IVy - 1 



for some constant ci > 0. If 5 > 1, then 



a c 



J(b + S)>-^= (1 y) 



which proves (2.22). A more careful computation for small 5 yields (2.23). 
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